Entanglement in the stabilizer formalism 
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Entanglement is an important, quantifiable physical 
resource of fundamental interest, rich with potential ap- 
plications in cryptography[l|, computation 2] , and con- 
densed matter systems Q. Bipartite pure state entangle- 
ment is the best understood 4] ; such states can be asymp- 
totically inter-converted at an exchange rate governed by 
the entropy of entanglement J3 of the original and trans- 
formed forms 0,13. A partial order on the space of such 
entangled states has also been found 8], leading to the 
discovery of ways to catalyze certain transformations us- 
ing other entangled states |9|]. 

However, because quantum states are generally impos- 
sible to describe concisely, e.g. an n qubit pure state 
may have 0(2"^) complex amplitudes, even well-defined 
measures such as the entropy of entanglement are hard 
to compute. 

The lack of efficiently computable entanglement mea- 
sures has also limited our understanding of the proper- 
ties of entangled quantum states shared between more 
than two parties. While the maximally entangled two- 
qubit singlet state plays the role of a "gold standard" 
for bipartite states, the three-qubit GHZ fails in this role 
for tripartite states [li|. In general, it is not known how 
to properly "price" multipartite entanglement, so inter- 
conversion relations are not well understoodlll|.ll2J. 

A wide class of interesting entangled states is the set 
used in quantum error correct ion |13j. cluster-model |l4j 
and fault-tolerant quantum computationjl^, and many 
cryptographic protocols such as secret sharing [l^. These 
are stabilizer states (also sometimes called graph states)^ 
and the study of their entanglement was introduced by 
Hein et a/.[l3], using a method of graphs. It was discov- 
ered that multipartite entangled states fall into a variety 
of equivalence classes, but the entanglement had to be 
quantified using the Schmidt measure 17], which is gen- 
erally computationally intractable. 

Here, we introduce a new method for computing the 
multipartite entanglement of any stabilizer state (includ- 
ing continuous variable stabilizer states such as coherent 
and squeezed states). Our measure of entanglement is de- 



fined for multipartite states, and is equal to the entropy of 
entanglement (up to a factor of two) for bipartite states. 
It can also be computed easily, requiring only a number 
of elementary operations which scales polynomially with 
the logarithm of the size of the Hilbert space. 

This result is made possible by the existence of efficient 
descriptions of stabilizer states, which require only 2n^ 
bits to specify an n qubit state |?/^). These numbers spec- 
ify the set of operators U in the Pauli group (i.e. tensor 
products of the identity / and pauli operators X, F, and 
Z) such that U\il)) = lip) (they stabilize |?/^)). These op- 
erators form a group S, generated by n operators, which 
we write as S = (^1,^2, •• • ^9n)' In terms of S^ we may 
express our main result in two parts as follows. 
Result 1: Entanglement of Stabilizer states: Just 
as the information content of a state \^|^)AB can be split 
into local information and correlations between A and B, 
the stabilizer S for \tp)AB can be split into a local sub- 
group Sa • Sb and a remaining subgroup Sab accounting 
for correlations. Sa {Sb) correspond to stabilizer opera- 
tors that act exclusively on A (5), as shown in Fig.^ For 
instance, the EPR state \0a Ob) + |1a 1b) is stabilized by 
S = {XX, ZZ), with Sab = S and Sa = Sb = {/}• Fur- 
thermore, the identification of these subgroups is simple, 
and require only 0{n^) steps for an n qubit state. 



SaI o.i 



Ib (1 < « < \Sa\) 






SbI Ia ® bj (^<j<\Sb\) 



FIG. 1: Canonical set of generators for a stabilizer group 
S{ip) with respect to a given partition {A, B} of the qubits. 
Sa and Sb contain the purely local information of \ip). Sab 
is generated by p pairs {gk^Qk) whose projections on A (or 
B) anticommute, but commute with all other generators of S 
including elements of other pairs. 



The first result of this letter is that the entropy of 
entanglement E{\il))) is given by 



E{m = -\SAB\. 



(1) 



where \Sab\ is the rank of Sab^ meaning the size of its 
minimal generating set. For the EPR state, \Sab\ = 2 
which correctly gives E = 1. For the three qubit GHZ 
state |000) + |111), where S = {XXX, ZZIJZZ), with 
respect to partition A = 12 and 5 = 3, we find Sa = 
{ZZI), Sb = {/}, and Sab = {XXX, IZZ), so we again 
obtain the correct result that E = 1, since l^^^l = 2. 

Importantly, this expression for E is easily computable; 
it requires only 0{ii?) operations. This is fundamentally 
because stabilizer states can be efficiently described in 
terms of the generators of their stabilizers. Furthermore, 
this stabilizer formalism give a constructive and efficient 
method to transform any bipartite stabilizer state by lo- 
cal unitary operations into E independent EPR pairs. 
These properties are proven below. 

Result 2: Multipartite entanglement: The stabi- 
lizer methods also apply to characterize the multipar- 
tite entanglement of stabilizer states. For an n qubit 
state l?/^) stabilized by S and split into k partitions 
A = {Ai, ^2, . . . , A/c}, we introduce a new, simple, mea- 
sure for multipartite entanglement. 



e^(IV^)) = n- 



n^. 



(2) 



where Sj contains the local operations of S that act as 
identity on the partition Aj. 

e^ is a true measure of multipartite entanglement. 
It is an entanglement monotone, meaning that it does 
not increase under relevant local operations and classical 
communications. For finer partitions, this entanglement 
measure is smaller, and it reduces to twice the entropy 
of entanglement for bipartite states. 

Finally in contrast to previously studied measures, it 
can be computed in 0{k-n^) steps. For example, Hein et 
al. utilize the Schmidt measure to characterize the en- 
tanglement of a class of graph states 17]. The Schmidt 
measure requires a difficult optimization for its computa- 
tion, limiting current studies to a small number of qubits. 
Graph states are also stabilizer states, and with this new 
stabilizer method, prior graph state equivalence classes 
with respect to local unitaries can be retrieved, using 
only simple manual computations; an example is shown 
in Figure El 

We now prove the above two parts of our result. 
Proof 1: The entropy of entanglement of a bipartite 



state 



lAB 



is defined as: 



E{\lP)) = -TY{pB\ogpB), 



(3) 




ei/2/3/4/5 = ^ 
e 12/345 = 2 
615/234 = 2 
^23/145 = 2 



^1/2/345 



= 1 



^12/3/4/5 — 1 




^1/2/3/4/5 = 
^'12/345 = 2 
^15/234 = 2 
^23/145 = ^ 
^1/2/345 = 1 
^'12/3/4/5 = 1 




^1/2/3/4/5 — 



f'12/345 - 

^15/234 

^23/145 

«l/2/345 

^1 2/3/4/5 = 



= 4 
= 4 



= 1 




^1/2/3/4/5 — 



^12/345 



= 4 



615/234 = ^ 
623/145 = ^ 
61/2/345 = 
612/3/4/5 = 



FIG. 2: Application of the measure ca to the classification 
of graph (stabilizer) states. The measure is shown only for 
the relevant partitions. A complete study would have to con- 
sider all partitions and relabelling of the qubits. From state 
(a) (GHZ) to state (d) (cluster), entanglement becomes more 
"localized" and robust against measurement of local opera- 
tors. 



where pB = Tr^ (|?/;) (?/;|). Since p = |'0)('0| is a stabilizer 
state, such that gp = p for all g G S, we may write it as 



p = ^Yl 



9' 



(4) 



g^s 



The partial trace over A thus gives the reduced density 
matrix 



PB 






2nB 



(5) 



Q^Sb 



where Sb is the subset of elements in S which are nonzero 
when traced over A. The entropy of pB is thus 



E{\^))=nB-\SB\. 



(6) 



SO the entanglement now reduces to computing the rank 
of 5b. 

This is accomplished most conveniently by using 
knowledge about the structure of the stabilizer S for 
\i^) AB with respect to the partition {A, 5}. Let Pa be 
the map that takes gA^QB ^ S onto gA^^B] this a pro- 
jection operator, such that Sb = KerP^. We construct 
a list of generators for S, by first including \Sa\ gener- 
ators ai Ib of Sa and | Sb \ generators I a ^ bj of Sb • 
Together, these generate the subgroup Sioc = Sa • Sb 
that we call the local subgroup of S. Assuming the list is 
completed with cab = n — \Sa\ — \Sb\ operators gener- 
ating the subgroup Sab^ then, S can be decomposed as 
S = Sa-Sb-Sab- Note that \Pa{S)\ =n-\SB\- Pa{S) 
is in general a non-abelian subgroup of the Pauli group. 

A good choice of operators to generate Sab can be 
found by studying the structure of subgroups such as 
Pa{S). First, it is helpful to define, for an arbitrary 



group G of Pauli operators, the compatibility index c{G) 
of G as the maximum rank of an abehan subgroup of G. 
Note that 1 < c{G) < \G\. For later convenience, also let 
the incompatibility index of G be p{G) = \G\ —c{G). The 
key insight into the subgroup structure (as illustrated by 
Fig. ^ is given by the following theorem: 

Theorem 1: The generators for stabilizer 5 of a bi- 
partite state can always be brought into the canonical 
form: 



S = {ai ^IbJa ®hj,gk,gk) ■ 



(J) 



where the first two subsets generate Sa and Sb-, and the 
last two generate Sab- These generators of Sab collect 
into p = p{Sab) anti-commuting pairs {gk^9k)^ where 
PA^dk) commute with all canonical generators of S ex- 
cept Qk^ and PA{gk) commute with all canonical genera- 
tors of S except gk . □ 

This theorem implies that a stabilizer state can be 
transformed into p independent Bell pairs by local uni- 
taries. A corollary of this is the relation \Sa\-\-'^P-\-\Sb\ = 
n. Since |5'a| +p < ^a and |5'b| +p < ^b, it follows that 



p = riA- \Sa\ = riB - \Sb\ = 



\Sab\ 



(8) 

\SbI 



It is also useful to know that since \Pa{S)\ = n 
p= \PA{S)\-nA = \PB{S)\-nB. 

Returing to our computation of the entanglement 
E{\ip))^ we now employ Eq.Q in Eq.Q and find that 
^(1^)) =p= \Sab\/2, as claimed in Eq.(P). 

The formalism used above shows that the problem 
of computing E{\ip)) reduces to the search of anti- 
commuting pairs in the projections on A or B of the 
generators of S. This takes 0{n'^ • min(n^,n5)) compu- 
tation time and uses 2n^ storage bits. Equivalently, one 
can compute the rank of Pa{S)^ which is the rank of a 
n X 2n matrix with elements in Z2 13j. 

The quantity E{\tlj)) has a particularly simple inter- 
pretation for graph states. In this case the group S has 
generators gj = Xj Yij. ^/c, where the product is over all 
nearest neighbors of j. An arbitrary element g ^ S can 
be written as ^ = YYj=i{9j)^^ ^^^ some binary vector 
X = (xi,...,Xn). Element g belongs to Sa under cer- 
tain circumstances. First of all we must have Xj — 
for all j G B. Then g acts as an identity operation on 
qubit j e B iS j has even number of neighbors k e A 
with x/c = 1. This requirement is equivalent to the Z2- 
linear constraint X]/cga ^j,k Xk = with F being the ad- 
jacency matrix of the graph. Thus elements of Sa are 
in one-to-one correspondence with zero vectors of an ad- 
jacency submatrix Tb,a between B and A. This proves 
that E{\tlj)) may also be computed as the binary rank of 
rB,A, for graph states. Bipartite entanglement in stabi- 
lizer states can thus be visualized as arising from graph 
edges crossing between the two partitions. Note, how- 
ever, that the number of such edges does not directly 



give the entanglement, unless the graph is first put into 
the canonical form given by Theorem 1. 

We now return to prove Theorem 1, in three steps. 

Lemma 1: c{G) > ^ {p{G) < ^) 

Proof: Denote by C a maximum abelian subgroup of 
G, generated by c{G) elements Cj, and by C the sub- 
group of G generated by \G\ — c{G) elements Ck such 
that G — {cj^Ck)' Up to a multiplication by elements of 
C, each operator Ck can be made to commute with all 
but one of the Cj. But then by the pigeon hole princi- 
ple, if c{g) < y , we can find ki 7^ /c2 such that Ck^ and 
c/c2 anti-commute with the same Cj, so that the product 
Cfei • c/c2 would commute with Cj (as well as with all the 
other generators of C), and hence C would not be the 
maximum abelian subgroup of G. D 

Lemma 2: We can choose the generators of G to be 
te}i..c(G) U {gj}i..\G\-c{G), where gj commute with ah 
operators except ^j, and gj commute with all operators 
except gj. 

Proof: From Lemma 1, we know that the generators 
can be organized as c{G) operators gj generating C{G) 
and \G\ — c{G) operators gk generating C{G)^ and also 
that each gj can be made to anti-commute with gj only. 
We now add a slight modification, recursively. Suppose 
Lemma 1 is obeyed if we keep only the first m anti- 
commuting pairs {gk-tQk)' Note that gm commute with 
all generators of G except gm^ and gm commutes with all 
gk^m and gk<m- If Qm+i and g^i do not commute, we 
redefine gm+i to be gm • ^m+i, so that up to this change, 
the Lemma is obeyed for the first tti + 1 pairs. D 

Note that the unpaired operators {^j, \G\— c{G) + 1 < 
j < c(G)} generate the center of G, the subgroup Z{G) 
that commutes with all elements of G (see Fig. ^. 



Z(G) 



C(G) 




C(G) 



FIG. 3: Structure of a subgroup G of the Pauli group. C{G) 
is a maximum abelian subgroup of S. It always includes the 
center Z{G) of G. 

The technical result that unravels the appropriate 
structure of S is the following: 

Lemma 3: The center of Pa{Sab) is trivial. 

Proof: Let z denote the rank of the center of Pa{Sab), 
and p its incompatibility index. We re-organize the 
generators of Sab so that their projection on A obeys 
Lemma 2. Note that their projections on B have a corre- 



spending structure, because the generators of Sab com- 
mute. Then taking into account the generators of Sa 
and Sb , we find \Sa\ ^ z ^ p independent commuting 
operators on A, and l^^l + 2: + p such operators on B. 
Therefore the fohowing inequahties must hold: 



\Sa\+ z+p <nA 


(9) 


\Sb\+z + p <nB 


(10) 


Sa\ + \Sb\+2z + 2p <n. 


(11) 



However, a simple generator count yields \Sa\ + \Sb\ + 
2p + 2: = n, so that z = 0. The construction used in the 
proof of Lemma 3 brings S into the canonical form of 
Theorem 1. D 

Proof 2: We now turn to the more difficult problem of 
finding a multi-partite entanglement measure for stabi- 
lizer states. This will be done by a generalization of the 
local subspace of S to the case of multi-partitions. 

Consider a /c-partition A = {Ai, ..., A/e} of the n 
qubits. We denote the projection on Aj by Pj for short. 
We define the subgroups Sj of S as Sj = {g e S, Pj (g) = 
/}, that is Sj is the kernel of Pj. We further define the 
local subgroup Sioc of S as 



^loc — 



n^. 



(12) 



In the bipartite case, Sioc = Sa' Sb- A qualitative differ- 
ence between the bipartite and multi-partite case is that 
the subgroups Sj might overlap. Therefore, it is harder 
to find a canonical structure for S when k > 3. Never- 
theless, the bi-partite case can be generalized to define 
a /c-partite entanglement measure e^ as in Eq.((2j, that 
is e^ = n — \Sioc\' For a bipartition, e^ reduces to cab 
which is twice the entropy of entanglement of Itp). 

To prove that e^ is actually an entanglement mono- 
tone, first note that each Sj and a fortiori Sioc are in- 
variant under local unitaries. Then note that the mea- 
surement of a local Pauli operator M can only increase 
\Sioc\' Simply, if M commutes with Sioc, then Sioc is 
contained in the new local subgroup. If not, then M re- 
places one generator of Sioc in the list of generators of 
the post-measurement state, but since M itself is local, 
the new local subgroup has not decreased in size. Finally, 
adding separable ancilla qubits to the system increases n 
and \Sioc\ by the same amount, and leaves the difference 
invariant. 

The entanglement measure e^ has also nice properties 
with respect to partitions. We say that a partition A is 
finer than a partition B {A ^ B) if every Ai is contained 
in a Bj. It is easy to see that if A ^ B^ then e^ < es- 
Simply, every Bj is a union of some A^, and therefore 
Sf^^ C Sj^^. Since each Sj can be found in O(n^) com- 



putational steps, the measure e^ can be computed in 
0{k • max(n^)) time, requiring 2n^ bits of storage. 

In summary, we have developed a mathematical for- 
malism to efficiently study the entanglement properties 
of stabilizer qubit states, which were already known to 
have an efficient classical description. Among other ap- 
plications, this formalism might be useful to study entan- 
glement in a quantum computation involving stabilizer 
codes. It could also be used to guide the construction 
of bipartite and multi-partite entanglement witnesses as 
combinations of stabilizer group generators. From a more 
fundamental point of view, it gives some insight into the 
problem of understanding multi-partite entanglement. 

As a final remark, we point out that this formalism 
is straightforward to generalize to qudits and continuous 
variable (CV) stabilizer states. For CV stabilizer states 
(including gaussian states), the Pauli group is replaced 
by the Heisenberg-Weyl group of displacement operators, 
but the formalism is the same. 

The authors gratefully acknowledge the hospitality of 
the Les Houches summer school, where this work was per- 
formed at the session on quantum information and quan- 
tum entanglement, directed by Jean-Michel Raimond 
and Daniel Esteve. This work was partially funded by 
the NSF Center for Bits and Atoms Contract No. CCR- 
0122419. 



[1] 
[2] 

[3] 
[4] 

[5] 
[6] 

[7] 
[8] 
[9 

[10: 

[11 

[12 
[13 

[14 

[15 

[16 

[17 



A. K. Ekert, Phys. Rev. Lett. 67(6), 661 (1991). 

A. Ekert and R. Jozsa, Proceedings of the Royal Society 

of London, Series A 356, 1769 (1998). 

G. Vidal et al., Phys. Rev. Lett. 90, 227902 (2003). 

V. Vedral and M. B. Plenio, Phys. Rev. A 57, 1619 

(1998). 

W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

C. H. Bennett et al., Phys. Rev. Lett. 76, 722 (1996). 

C. H. Bennett et al., Phys. Rev. A 53, 2046 (1996). 
M. A. Nielsen, Phys. Rev. Lett. 83, 436 (1999). 

D. Jonathan and M. B. Plenio, Phys. Rev. Lett. 83, 3566 
(1999). 

W. Diir, G. Vidal, and J. L Cirac, Phys. Rev. A 62, 
062314 (2000). 

C. Bennett et al., Phys. Rev. A 63, 012307 (2000). 

W. Diir and J. L Cirac, Phys. Rev. A 61, 042314 (2000). 

D. Gottesman, Ph.D. thesis, California Institute of Tech- 
nology, Pasadena, CA (1997). 
R. Raussendorf, D. E. Browne, and H. 
Rev. A 68, 022312 (2003). 
J. Preskill, Proc. Roy. Soc. A: Math. 
454, 385 (1998). 

R. Cleve, D. Gottesman, and H.-K. Lo, 
83, 648 (1999). 

M. Hein, J. Eisert, and H. J. Briegel, 
062311 (2004). 



J. Briegel, Phys. 

Phys. and Eng. 

Phys. Rev. Lett. 

Phys. Rev. A 69, 



